CERN-TH/98-143 



hep-ph/9805445 



Weyl— van-der-Waerden formalism 
for helicity amplitudes of massive particles 

Stefan Dittmaier 

Theory Division, CERN 
CH-1211 Geneva 23, Switzerland 



Abstract: 

The Weyl-van-der-Waerden spinor technique for calculating helicity amplitudes of 
massive and massless particles is presented in a form that is particularly well suited 
to a direct implementation in computer algebra. Moreover, we explain how to exploit 
discrete symmetries and how to avoid unphysical poles in amplitudes in practice. The 
efficiency of the formalism is demonstrated by giving explicit compact results for the 
hehcity amplitudes of the processes 77 f f, ff — >• 777, /x"/^"*" — > //7- 
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1 Introduction 



The calculation of scattering amplitudes and cross sections in lowest order of perturbation 
theory is one of the standard problems of elementary particle phenomenology in practice. 
Although this task is very simple for 1^2 particle decays and 2 — 2 particle reactions, 
the situation can become arbitrarily complicated for more particles in the final state, as 
the number iV^j of contributing Feynman diagrams rapidly increases. Squaring amplitudes 
and using completeness relations for the wave functions does not represent an appropri- 
ate approach if Nd is large, because one would get Nf^ contributions. A reasonable way 
to avoid an explosion of the algebraic work consists in using spinor formalisms for the 
evaluation of amplitudes for definite polarization configurations. The squaring of the am- 
plitudes and the spin summations are then performed numerically. Of course, the number 
of different polarization configurations also grows for more and more external particles, 
but in practice this fact is much less problematic than squaring amplitudes. In particular, 
the number of independent matrix elements can be reduced by exploiting discrete symme- 
tries. Although the commonly applied spinor formalisms rely on the four- dimensionality 
of space-time, they are also useful for the calculation of higher-order corrections if the 
singularities in dimensional regularization are split off and controlled separately. 

Various spinor techniques have been proposed and successfully applied in recent years. 
The first versions proposed a clever choice of photon and gluon polarization vectors, 
which was dictated by the momenta of the attached fermions. This choice forces a lot of 
terms in the calculation to vanish, and relatively compact amplitudes result. Although 
the initial restriction to massless particles could be overcome [0] , all those formalisms have 
in common that the usual Dirac algebra is still present. The final results for amplitudes 
are given in terms of standard products Ui{pi)u2{p2) of Dirac spinors. This form is useful 
for numerical evaluations, but too involved for further analytical purposes. 

It is more elegant and transparent to express all needed ingredients, i.e. momenta, 
polarization vectors, and Dirac spinors, in terms of a single mathematical object. This 
leads us to two-component Weyl-van-der-Waerden (WvdW) spinors [^, which form the 
fundamental representations of the Lorentz group, and the related spinor calculus (see also 
Ref. [0). Since WvdW spinors are closely related to light-like vectors, the application of 
the WvdW spinor calculus is extremely simple for massless particles. Such applications 
can be found, e.g. in Ref. |]^. The WvdW technique can also be used for massive particles 
if time-like momenta are decomposed into two light-like ones. This procedure was applied 
in Refs. [^, for the calculation of bremsstrahlung corrections to production processes 
of massive weak gauge bosons, leading to rather compact results. However, the WvdW 
spinor formalism for massive particles is not yet well documented in the literature. In 
connection with the actual applications presented in Refs. |^, the formalism was only 
briefly sketched, as far as was relevant for the special cases under investigation. More 
complete presentations of the technique for massive particles are either not suited for 
direct practical use |Q or appeared in unpublished studies |P, [10|. 

In addition to the more technical motivation above for considering amplitudes with 
definite polarization configurations, one should add some remarks on the physical role of 
polarized massive particles. In contrast with massless particles, the spin orientations of 
massive particles do not transform covariantly under Lorentz transformations; in partic- 
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ular, the property of being a helicity eigenstate is frame-dependent. Nevertheless, the 
polarization states of a massive particle also carry valuable physical information about 
the structure of its interaction. For instance, the longitudinal modes of the massive weak 
gauge bosons serve as a window into the scalar sector of the electroweak theory, and the 
reconstruction of the top-quark spin from its decay products reveals information about its 
static properties and form factors. For various aspects of physics with polarized massive 
particles, see Ref. [|TT| and references therein. 

The aim of this paper is threefold. First, we give an introduction into the WvdW 
spinor technique for massive and massless particles. The presentation is deliberately held 
at a basic and detailed level, in order to facilitate the practical use of the method. In 
particular, we present all ingredients needed for an implementation in computer algebra. 
As wave functions for massive spin-| and spin-1 particles, we take the helicity eigenstates 
formulated in Ref. |jTO|, which appear to be very simple. The wave functions and the 
momenta of the corresponding massive particles are composed of the same set of auxiliary 
spinors, leading to simplifications in the calculation of amplitudes. Secondly, we show 
how to exploit the advantages of the formalism in practice. We explain the use of discrete 
symmetries, in order to reduce the algebraic work, and give a prescription for avoiding 
unphysical poles in amplitudes. In the presented examples, we treat the polarizations of 
massive fermions in a generic way, i.e. explicit results for definite helicity configurations 
follow from a few compact generic amplitudes upon setting auxiliary spinors to specific 
values. Thirdly, we give explicit results for helicity amplitudes of phenomenologically 
relevant processes involving massive particles. We explicitly treat the reactions 77 ff, 
ff — >■ 777, //"/i"*" — >■ //7 in the framework of the electroweak Standard Model, and show 
how these results can be carried over to other processes such as 77 ff'J, e~7 £"77, 
e~7 —>■ e~e~e"'" by using crossing relations. 

The outline of the paper is as follows. In Sect. |^ we set our basic conventions for 
the WvdW spinor calculus and introduce helicity eigenstates for spin-| and spin-1 parti- 
cles. The necessary ingredients for calculating helicity amplitudes are described in Sect. |^, 
specifically containing a prescription for translating Feynman rules into the WvdW for- 
malism, a formulation of discrete symmetries for helicity amplitudes, and a prescription 
for avoiding unphysical poles in amplitudes. The physical applications are presented in 
Sect. |, and Sect. ^ contains a summary. 



2 Basic definitions 
2.1 Spinors 

The basic philosophy of the WvdW formalism is to reduce all mathematical objects that 
belong to higher-dimensional representations of the Lorentz group to the two-dimensional 
irreducible representations D(|,0) and D(0, |). According to these representations we 

distinguish covariant and contravariant WvdW spinors ipA and V'"^, respectively. Indices of 
WvdW spinors are denoted by capital letters throughout. The transition between the two- 
dimensional representations, which are non-equivalent, is achieved by complex conjugation 
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and a similarity transformation. Complex conjugation is consistently indicated by dotting 
(undotting) indices, i.e. 

V'i = {^a)\ = {^^T- (2.1) 

The similarity transformation is mediated by the antisymmetric 2x2 matrix e = io"^, 
where cr" (a = 1, 2, 3) are the standard Pauli matrices, 

e^^ = e^^ = e^^ = e^^=(_°^ . (2.2) 

The matrix e defines how to raise and lower spinor indices, 

V'^^e^^^B, ^^ = e^^VB, ^iPA = V^''eBA, ^A = i^''^BA- (2-3) 

Moreover, e provides the set of Clebsch-Gordan coefficients that project the product 
representations D(|,0) ® D{^,0) and i^(0, |) ® -D(0, |) onto the trivial representation. 
For this reason, it is possible to define a Lorentz-invariant spinor product through 

{(PiP) = = 01^2 - </>2V'l, ('/•V')* = = (0lV'2 - 02^l)*, (2.4) 

and e is called spinor metric. By definition, the spinor product is antisymmetric: 

(0^) = -m, (00) = 0. (2.5) 

The fact that an antisymmetric tensor built of two-dimensional objects cannot have a 
rank higher than 2 implies a Schouten identity for the spinor metric, 

^AB^CD ^ ^AC^DB ^ ^AD ^BC ^ (2.6) 

which in terms of spinor products reads 

{4>^){iv) + + = 0- (2.7) 

2.2 4- vectors 

Minkowski 4-vcctors belong to the representation D(|, |) = D{^,0) -D(0, |) of the 
Lorentz group. The transition of the usual form of a 4-vector k'^ = [k^, k) to the spinor 
representation D(|, |) is provided by the matrices 

a'^'^^^{a',a), a^^- ^ {a', -ct), (2.8) 

consisting of the 2-dimensional unit matrix (7° and the Pauh matrices cr". Each 4-vector 
k^ is related to a 2x2 matrix 

= (^i^ip A;O-A;0' ^^'^^ 

which is Hermitian if the components of k'^ are real. The rules for dotting, undotting, rais- 
ing, and lowering spinor indices also apply to the indices of the a matrices; in particular, 
we have 

<B = 'CA eDB, = {a^^^y. (2.10) 



We note that the coefficients of the transpose K"^ of a matrix K read K^^ if the ones of 
K are denoted by K^^, i.e. transposing a matrix interchanges the spinor indices without 
moving the position of the dot. Thus, the Hermiticity of the a matrices is expressed by 



a a 



(2.11) 



The a matrices obey the important relations 



^Ab ^f^,CD 



2k-p = k^2g^^p, = k^a^^^a'^^^^p, = Kj^^P^^, (2.13) 



(2.12) 

The first of these relations translates the Minkowski inner product of two 4-vectors k^ 
and p^ into 

and the second one implies 

K^j^K^^ = kH%. (2.14) 

In order to reduce terms involving a matrix Kj^^^ to spinor products, it is necessary 
to express Kj^^ in terms of spinors. For a real 4-vector the matrix Kj^^^ is Hermitian and 
can be decomposed into its eigenvectors rii^A {i = 1, 2) and eigenvalues Aj: 



K 



AB 



XiUiArii^B, Ai,2 = k^ ± |k| 



i=l,2 



e ^ cos ^ 

. a ^ 

sm ■ 



n2,A 



sm 

-e+''^ COS; 



:2.15) 



where Q and denote the polar and azimuthal angles of k = |k|e, respectively, 

/ cos sin Q \ 
e = sin sin Q 
\ cos 9 J 



(2.16) 



For time-like vectors (/c^ > 0) it is often convenient to include the eigenvalues Aj in the 
normalization of the eigenvectors, resulting in 



^Ab — X! ^j,A '^i,B, 



l^i,A 



A,: Hi A- 



[2.17) 



«=1,2 



The phases of the rii^A are chosen such that the orthonormality relations read 

(njrii) = 0, (risrii) = -(riiris) = +1. (2.18) 
Note also the relations for the eigenvalues Xf. 

det{K^j,) = A1A2 = (/t2/€i) = v/V^ = v^- (2.19) 

The special case of a light-like vector (/c^ = 0) is of particular importance. In this case 
the eigenvalue A2 of ( |2.15| ) vanishes, and the matrix Kj^^ factorizes into a single product 
of two spinors 



^Ab — ^^^ijl,AB — ^A^B, 



'2k^ ( ^ . 2 
sm2 



(2.20) 



In this context, kA is called a momentum spinor. 

Finally, we remark that the decomposition (|2.15| ) is a very convenient, but not unique, 
possibility to express a 4-vector with fc^ 7^ in terms of WvdW spinors. Any splitting 
of k^ into two light-like 4-vectors yields a decomposition into spinors, since light-like 
vectors factorize, as seen above. For instance, choosing an arbitrary light-like 4-vector 
[a? = 0) with a ■ k ^ 0, and defining 



1.2 

a = -, b'' = k>'- aaf", (2.21) 

2a -k' ' ^ ^ 



yields a possible decomposition k^^ = aa^ + h^. In terms of WvdW spinors this corresponds 
to an arbitrarily chosen spinor a a with K(jj^a^a^ 7^ 0, leading to the decomposition 

K ■ P 
Kj^j^ = a aj^ttB + hjpB with 6^ = , . a = — p— . (2.22) 

3 Wave functions for helicity eigenstates 
3.1 Spin-| particles 

Dirac spinors \1/ belong to the representation -D(|, 0)©D(0, |) of the Lorentz group. Thus, 
in terms of WvdW spinors, they are represented by 

(3.1) 

The WvdW formalism consistently employs the chiral representation of the Dirac matrices 

■^'-[J.^s t). 7'^i7V-,V^(;° _°„). (3.2) 

We are interested in plane-wave solutions \& = exp{=piA;x}\l/[^^'' of Dirac's equation 

(i^-m)^ = 0, (3.3) 

which describe the propagation of free Dirac fermions and antifermions, respectively. In- 
serting ( |3.1| ) and ( p.2| ) into Dirac's equation, we obtain the coupled pair of Weyl equations 

ir^^#)'^ = ±m0g, ir^>^5 = ±m#)'^ e=m\ (3.4) 



Upon applying the decomposition ( |2.17|) of K^^, the following solutions can be easily 
constructed: 

The corresponding adjoint spinors \E' = \I^^7o read 

^i? = '^i.a) , = {f^t ±^^2^) ■ (3-6) 
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The solutions obey the standard normahzation 

= 2A:° 5.,, ^i?vl/W = ±2m5.±5,,. (3.7) 

Moreover, it is straightforward to check that they indeed form eigenstates of the hehcity 
projector 

S^ = Kl±7V.), -^ = ^|j-^^°f|- (3-8) 

More precisely, projects onto "^^f^^ and \1'^ 2^ o^^^o ^^'^ ^1^2*' i-^- ^i'l'' describes 
a right-handed fermion, \E'[. 2^ a right-handed antifermion, etc. 



For massless fermions the Weyl equations (p^ ) decouple, and the helicity eigenstates 
take the simple form 



7(±) /n N 7f7(±) 



^M=(0,fc^), v&S = (A:^0). (3.9) 



For the decomposition (|2.22| ) of K^^, also simple plane- wave solutions exist, namely 



which are, however, not related to definite helicity states. 
3.2 Massive spin-1 particles 

Spin-1 fields V'^ transform like ordinary 4-vectors under Lorentz transformations. This 
means that polarization vectors e'^ for such fields can be related to 2x2 matrices ej^^ in 
the same way as described in ( |2.9| ) for a general 4-vector fc^. However, one should keep in 
mind that polarization vectors need not be real so that e^^ is not Hermitian in general. 
The free field V'^ of a massive vector particle is governed by Proca 's equation 

{d^ + m')g^"' - d^d"] V, = 0, (3.11) 

which for m 7^ is equivalent to the Klein-Gordon equation, in conjunction with the 
transversality condition: 

{d'^ + m^)V^ = 0, d-V = 0. (3.12) 

With the ansatz = exp{'^ikx}e^{k) for particles and antiparticles, respectively, we 
arrive at 

fc2 = k^e^{k) = iK^^e^sik) = 0. (3.13) 

For m 7^ there are three linearly independent, space-like polarization vectors e^{k), 
which are usually orthonormalized according to 

enk)sl,ik) = \e,^AB{k)ef''{k) = z,j = 0, ±, (3.14) 
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where e*"^^ = e*a'^'^^ . A helicity basis for the s^{k) is given by 



em 



V2 



(0, — COS 6 cos ± i sin 0, — cos 6 sin =i= i cos 0, sin ( 



— — , cos sin 6', sin sin 6*, cos 9 



(3.15) 



where is related to as given in (3^). In terms of WvdW spinors this spin basis reads 



B-, 



'0,AB 



{k) 



m 



(3.16) 



where the spinors n^^^ and defined in (|2.15|) and (|2.17|) , respectively. Note that care 

has to be taken when dealing with conjugate polarization vectors, which, in particular, 
occur for outgoing spin-1 particles in transition amplitudes. In such cases the transition 
from the 4- vector e*'^ to the 2x2 matrix e*j^^ upon contraction with the a matrices has to 
be consistently performed for the conjugated polarization vector, i.e. e*j. 



AB 



e ,a' 



IJ- AB 



which 

in general is different from taking the complex conjugate of sj^^. For clarity, we give also 
the helicity basis for an outgoing spin-1 particle: 

1 . 



,AB^ 



m 



(3.17) 



which correspond to the conjugate polarization vectors e* of given in (|3.15|) . From 
( p.l6| ) and (|3.17|) we obtain the relations 



'i.AB 



{k) 



,BA 



{k) 



-i,AB 



{k) 



^i,BA 



{k). 



(3.18) 



Finally, we give a simple spin basis for the decomposition ( |2.22| ) of -ft'^^. 



-+,ABiP)- ^^^^y > --,AB 

which is, however, not related to helicity eigenstates. 



e 



^o,ab(p) = —{b^bB-a a^as), (3.19) 



3.3 Massless spin-1 particles 

In contrast with the case of spin-i particles, the zero- mass limit for vector bosons confronts 
us with a new physical situation. For m — *• the polarization vector does not exist, 
since the longitudinal polarization does not yield a physical state for massless spin-1 
particles. The corresponding degree of freedom manifests itself in the arbitrariness of 
gauge for the two polarization vectors e±. In the language of WvdW spinors this fact is 
expressed very elegantly: 

_ ^9+,A kB V2k^g_^B 



{g^ky ' ' (g^k) 
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Outgoing fields: 


— 1-1? 




— 


• 4 










Propagators: 


B k A i(^ + m/) ^ i fmf6^ 

k'^-mj k'^ - mj \ K"^^ ^f^i J 




AB CD -igfiu -2ltj^(jeBD 



Table 1: Prescriptions for translating ordinary Feynman rules for external fields and 
propagators into the WvdW formalism. 

where g±^A denote arbitrary spinors with {g±k) ^ 0; they are called gauge spinors. The 
difference of two matrices ^iAsi^) differently chosen gauge spinors is proportional to 
the momentum matrix kj^ks, as can be easily checked by applying the Schouten identity 
( [^.61 ). The transversality condition k ■ e± = is still fulfilled for any g±^A, i-e. the Lorentz 
gauge is maintained for the e± of ( p.20|) . Transition amplitudes do not depend on the 



choice of gauge spinors. This means that the g±^A can be appropriately chosen to simplify 
the algebraic expression. Note that it is even possible to take different sets of gauge 
spinors for each helicity channel. 



4 Calculation of helicity amplitudes 
4.1 Feynman rules 

Feynman rules are usually written down in terms of 4-vectors and Dirac matrices. Here 
we describe how such Feynman rules can be directly translated into the language of 
WvdW spinors. Obviously this procedure is much more practical than going back to the 
underlying Lagrangian and introducing WvdW spinors there. 

Given any Feynman graph, the trick is to contract each vector-boson leg of a vertex 
with the identity = ^a^^a'^^ and to shift the factor cx^^ to the vector-boson propagator 
or to the external wave function of the vector boson that is attached to this vertex. 



Moreover, for each occurring Dirac matrix the chiral representation (|3.2| ) has to be used 



In particular, the unit matrix 1 in the Dirac space and a slashed quantity read 

(4.1) 





^AB 
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A 71/2 



le 




\eCv1V2V3 



AB 

CD 




ieCsviV^g^" 



AC^BD 



^eCsviV2^ e 



AB . ; 

ieCy5,&(/^i-^2)'' ^ ^eCvs,S2{Ki-K2)^'' 



AS 

CD 




EF 



i ^2 



Viy2V3y4 



2^AC^BD^EG^FH 



^AE^BF^CG^DH 



_ ^AG^BH^CE^DF 



L'ViV2S3S4^ e 



Table 2: Prescriptions for translating ordinary Feynman rules for vertices into the WvdW 
formalism. Scalar fields, fermion fields, and vector fields are generically denoted by 5", /, 
and V, respectively, and u!± — |(1 ± 7^) are the chirality projectors. 
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respectively. Tables |l] and |^ summarize the necessary changes of the generic Feynman 
rules for the electroweak Standard Model in the 't Hooft-Feynman gauge. The explicit 
values for the couplings can, for instance, be found in Ref. |]12|. The new Feynman rules 
for other conventions or different models can easily be worked out by the reader. 

Once the Feynman rules are settled, it is very easy to explicitly write down helicity am- 
plitudes for any process at tree level. Expressing all particle momenta in terms of spinors, 
as described in Sect. each helicity amplitude is reduced to an algebraic expression 
in terms of antisymmetric spinor products {(p'lp) after all spinor indices are contracted. 
These spinor contractions can be performed like usual contractions of Lorentz indices, 
apart from taking care of the antisymmetry. The form of the amplitudes obtained this 
way is already well suited to numerical evaluations. 

4.2 Discrete symmetries 

Discrete symmetries relate helicity amplitudes of different processes or of one and the 
same process. Thus, they either provide convenient cross checks for results or they allow 
for a reduction of the algebraic work by generating various amplitudes from a generic set 
of amplitudes. In the following we show how to derive the relations implied by crossing 
symmetry, parity, and CP symmetry, which are the most important discrete symmetries 
in practice. 

(i) Crossing symmetry: 

Crossing symmetry transforms an incoming particle into the corresponding outgoing 
antiparticle, or vice versa. Denoting the momentum of a given particle by the 
inversion k ^ —k can be consistently obtained by substituting 

^i,A ^ ~^i,Ai ^i,A ^ ~\~f^i,A (4-2) 

in the decomposition (|2.17|) of /c, i.e. by inverting the contravariant parts only. In or- 
der to relate helicity amplitudes, it is also necessary to consider the relation between 
the wave functions of the respective incoming and outgoing fields. Inspecting the 
explicit form of the wave functions for the helicity eigenstates of Sect. |^, one finds 
that the substitution ( |4.2| ) transforms incoming (outgoing) particles into outgoing 
(incoming) antiparticles with reversed helicity modulo sign change. Specifically, 
there is a global factor —1 for each crossed spin-1 field and a factor ±sgn(cr) for 
each incoming/outgoing fermion with helicity a that is involved in the crossing. 

(ii) Parity: 

If parity is a symmetry, every helicity amplitude, up to a phase factor, agrees with 
the corresponding amplitude with opposite helicities, after the spatial parts of all 
momenta are inverted. At the level of WvdW spinors, the inversion of the spatial 
parts of momenta is connected with the interchange of the spinors ipA and ip"^. It 
turns out that each helicity amplitude, up to a global sign factor, agrees with the 
complex conjugate amplitude with opposite helicity configuration. More precisely, 
if iA^((Tj; \j) is an amplitude involving n fermions and n antifermions with helicities 
(Ti, and ny vector bosons with helicities Aj, we get the relations 

M{-a,--\j) = {-lfsgn{ai---an+n)M{ai-\jY (4.3) 
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for tree-level amplitudes. If parity is violated, these relations for amplitudes remain 
valid if left- and right-handed couplings are appropriately substituted. The explicit 
derivation of ( [4.3| ) and the modifications for broken parity can be found in the 
appendix. 

(iii) CP symmetry: 

CP is the product of parity, which is explained above, and charge conjugation, which 
interchanges particles with the respective antiparticles. Therefore, CP symmetry 
leads to relations between the helicity amplitudes of a given process and the complex 
conjugate helicity amplitudes with reversed helicities of the process involving the 
respective antiparticles. We consider an amplitude ij^^af, Xj) and the one for the 
CP-related process, iAi {ai] where the respective helicities for the fermions and 
vector bosons are assigned as above. If CP is an exact symmetry these two matrix 
elements are related by 

M-a,;-A,) = (-l)"^+"^"+*^^"sgn(ai---a„+^)A<(a,;A,)* (4.4) 

at tree level, where riin + ^in is the total number of incoming fermions and an- 
tifermions. The explicit derivation of these relations and their modifications if CP 
is violated are given in the appendix. 

The actual use of these symmetries will be illustrated when we consider helicity amplitudes 
of concrete processes, in Sect. |[ 



4.3 Remarks on the choice of gauge spinors 

Since a helicity amplitude does not depend on the actual insertions for gauge spinors, they 
can, in principle, be chosen arbitrarily. Usually a gauge spinor is chosen such that as 
many terms as possible in an amplitude vanish. However, as can be seen from ( ^.20 ), an 



arbitrary choice of gA in general leads to unphysical poles at the zeros of (gk), in individual 
terms that contribute to an amplitude. Of course, the unphysical pole drops out in the 
final result, but this cancellation can cause instabilities in a numerical evaluation. 

Unphysical poles can, for instance, consistently be avoided by setting gA = ^2,A5 where 
n2,A is related to the momentum k of the massless spin-1 particle as specified in ( p^.l5| ). 
This choice is identical with ( |3.16|) for the transverse modes of a massive spin-1 particle. 



The drawback of this choice is that no algebraic simplifications result. 

Another possibility to avoid numerical problems that are due to unphysical poles is to 
cancel such poles analytically, before the numerical evaluation. In general this task can 
be extremely cumbersome, and the additional work devalues preceding simplifications. 
In the following we give a very convenient choice of gA, yielding desirable simplifications 
without leaving uncancelled unphysical poles. 

Consider a process with an incoming photon of momentum k and an outgoing fermion 
/ of momentum p; the momentum matrix for p is denoted by 

Pab = E f^i,A'^i,B- (4.5) 

i=l,2 
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Figure 1: Subdiagrams containing the fermion propagator denominator [{p — k)"^ — m^j^^ 
for processes with an incoming photon of momentum k {k"^ = 0) and an outgoing fermion 
or antifermion of momentum p {p^ = mj). 



We choose the same gauge spinor 



(4.6) 



i=l,2 



for both photon hehcities. The problematic denominators contained in the photon polar- 
ization vectors ( p.20|) are given by 



k^Pj^j^k^ 



This has the form of the inverse propagator {p — kY 



2{p-k). (4.7) 

roj = —2{p ■ k) that appears if the 
incoming photon directly couples to the outgoing fermion / so that no unphysical pole 
is introduced at all. The subdiagram containing this propagator is shown in Fig. |I]a. In 
Feynman graphs that contain this subdiagram, we obtain the factor [p ■ fc)~^, one power 
is due to the propagator, another is due to the polarization vector. For small [p ■ k) we 
again get an unwanted cancellation of one power in {p- k). However, this cancellation can 
be easily performed analytically. Denoting the Dirac spinor of the outgoing fermion with 
helicity a, generically, by 



^17 = (^^,0c) with (0,^/^) 



(ki,-K2) for i = 1 (a = +), 
{k2, Ki) for z = 2 (o- = -), 



(4.^ 



and using the wave function e^^ for the photon with helicity A, the subdiagram reads 



Qfe (^,c ^\( ^Pc\^ ( ^f^D {p-k)DE 



2{p ■ k) 



(V^"", 0c) 



AB 



(p-k) 



DE 



(4.9) 



where Qj is the relative electromagnetic charge of the fermion /. Upon inserting ( [4.5| 
([4.6|), and the exphcit expressions ( |3.20| ) for e^^^, we get 

Qfe 



{kip) [2{p ■ k)k'' + ruf (0V^)A;^ - + ,0 . (4.10) 



2V2{p ■ ky 

Here we have used that P^^ = 4>j^(f)B + and {(pip) = fo^' both fermion helicities. 



The expression within square brackets in ( [4.10| ) vanishes according to Schouten's identity 



12 



(2.6) so that one factor (p ■ k) cancels. Performing a similar calculation for Tl-^ (a), we 



find the simple results 

Considering an outgoing antifermion instead of a fermion, as shown in Fig. |l]b, we get 



where the Dirac spinor for the antifermion is generically denoted by 

with (^.VO^I*"""^' for. = l(a = -). 



From the results ( [4.11j ) and ( 4.12 ) for the subdiagrams (a) and (b) of Figs. |l], respectively, 



all similar subdiagrams involving outgoing photons and/or incoming (anti)fermions follow 
by crossing, as described in the previous section. 

The above procedure also works if the (anti)fermions in Fig. |I] are replaced by charged 
gauge bosons, e.g. by the W boson of the Electroweak Standard Model. In this case, how- 
ever, one has to take into account the contribution of the associated would-be Goldstone 
boson, which appears on internal lines (vertical lines in Fig. |I|). 

5 Applications 

In order to illustrate the actual use of the formalism described above, we calculate full sets 
of helicity amplitudes for some processes with massive particles and photons in lowest or- 
der. The presented results have been derived analytically and simplified as far as possible; 
this demonstrates the analytical power of the technique. In addition, the amplitudes have 
been evaluated by performing the spinor contractions automatically in Mathematica |l^ , 



which is very simple. The results of both approaches have been compared numerically. 

The following examples only include massive external fermions. For an application 
of the helicity basis for massive spin-1 particles we refer to Ref. J^, where the radiative 
processes e~7 W~z/e7,e~Z7 are discussed. 

5.1 The process 77 — > ff 

For illustration we start with the simple QED process 

7(^1, Ai) + ^{h, A2) fip, a) + fip', a'). (5.1) 

The momentum matrix and the wave function for the outgoing fermion / are denoted as 
in (|4.5|) and (O), respectively. For the outgoing antifermion we define 



1=1,2 
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(a) (b) 

Figure 2: Born diagrams for 77 — //. 

and 

Since we want to make use of the results of the previous section, we set the gauge spinor 
gn,A for the nth photon to 

gn,A = PaB^h = XI l^iM'^i^n)* , n= 1,2. (5.4) 

i=l,2 

In lowest order the two Feynman graphs of Fig. ^ contribute. Let us first calculate the 
contribution A^(a) of diagram (a) to the full helicity amplitude We identify the photon 
7(/c, A) in the subgraph shown in Fig. |l|a with 7(^1, Ai) and denote the contribution of 
this subgraph by Tl^J {a). M.{a) then follows upon multiplication of T^^^{a), which is 
given in ([4.11|) , with the second 7// vertex and the wave functions for the second photon 
and the antifermion, 

i^(.)(Ai, A2,a,a') = r,^y(a)(-iQ,e) (^^J^^^ '^^^^^ ) (^'^) ^^'^^ 

Note that we follow the common practice to split off a factor +i from helicity amplitudes 
throughout. It is trivial to carry out the spinor contractions in (|5.5|) for the explicit 
insertions of T^^J [a) and £2GH- '^^e contribution Ai{h) of the diagram in Fig. ^ follows 
from M.(a) by consistently interchanging the incoming photons. This leads us to the final 
result for the full helicity amplitudes 7V1 = (a) + (b) : 

2{p ■ ki){p ■ k2) ^ ' 

^(-,+,a,a') = ^(+,-,a,a')^ . (5.6) 

fcl^fc2 

Here we have introduced the abbreviation 

{hPk^) = ki^^P^''kn,B = E {kl^^^y{k^K,), l,n = 1,2. (5.7) 

i=l,2 
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Moreover, we mention that Schouten's identity ( |2.6| ) has been used in order to get compact 
results for Ai = A2 = ±. 

As explained in Sect. parity relates amplitudes with opposite-helicity configu- 
rations through complex conjugation, modulo a sign factor. According to (^4.31 ), these 
relations explicitly read 

Mi-Xu-X2,-cr,-a') = -sgn((Ta')A^(Ai,A2,a,cr')*, (5.8) 

consistent with the above results. Alternatively, these relations can be used to generate all 
hehcity amplitudes from the generic results for A^(+, +, a, a') and A4{+, —,a,a') given 
in ( ^.6|) . The process 77 — > // is also CP-symmetric. Taking into account that the CP 
transformation interchanges / and /, Eq. ( [4. 41) implies 



M{-Xu -A2, -a', -a) = sgn{aa')M{Xu A2, a, a')* . (5.9) 



where the substitution p p' also includes the interchange of the respective spinors Hi 
and k'. 



5.2 The process // — > 777 and related reactions 

Next we consider fermion-antifermion annihilation into three photons: 

f{p,(r) + f{p',a') 7(^1, Ai)+ 7(^:2, A2)+ 7(^3, A3) 



(5.10) 



The momentum matrices P and P' for the momenta p and p', respectively, are again 
decomposed into the respective spinors A and k^^, as defined in ( |4.5|) and (|5.2|) . The 
Dirac spinors are generically assigned by 



p 



p' 



with the actual insertions 



{1^2, 



for cr = +, 
for cr = — , 



(</>', ^') 



(^'^0i) 



(K^, K2) 



(5.11) 



for a' 
for a' 



(5.12) 



For all helicity configurations we have {(pip) = (V'V) = ""^z- The polarization spinors 
A^{ki) for the outgoing photons are defined as in ( 3.20 ). Following the strategy of 
Sect. ^]3| and "crossing" the results ([4.11|) or (|4.12| ) for the subdiagrams of Fig. Ill, the 
actual calculation of the six diagrams for // '-yyj (see Fig. ^) is rather simple. In order 
to minimize the number of generic amplitudes that have to be calculated, we first give 
the relations that follow from discrete symmetries. Parity, CP, and Bose symmetry imply 

M{-a, -a', -Ai, -A2, -A3) = -sgn(aa')A^(a, a', Ai, A2, A3)*, 

M{-(t', -a, -Ai, -A2, -A3) = -iign{aa')M{a, a', Ai, A2, A3 

M{(y, a\ Ail, Ai2, Aig) = M{a, a\ Ai, A2, A3 



(5.13) 
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■fV\/\/V 



irv\/W 



+ five permutations of tlie pliotons 



Figure 3: Born diagrams for // 777- 



Therefore, for any given {a, a') combination it suffices to calculate the helicity amplitudes 
with (Ai, A2, A3) = (+++), (+H — ), from which all other amplitudes follow by ( |5.13|) . 
These generic results can be expressed in a very compact form: 



M{a, cr', Ai, A2, A3) 



(5.14) 



where 

A, 



[a, a 



mf— —— ikiP' Pki) + cyclic permutations in /ci, A;2 5 ^3, 

[p' ■ ki) 



yl++_(a,a') = nij f^^^'^ {{k,P' Pk,)* {H') + 2(p ■ k,){k,cl^){k,^') 
{k2Pk^) r 



{p' ■ ki 



{k,p'Pk,){{k2tpr{k,tp') + {k2<j)r{h<p)) 

+ 2{p-ki){ks(j)){k,<f)')*{k,k2y] + (fci^fcs). 



(5.15) 



In addition to the shorthand {kiPkn) of (ISTTI) we have introduced the useful abbreviations 
{hP'Ph) = k,^j,P'^^PcBk? = (klK[Y{K,K[){K,h)\ / = 1,2,3. (5.16) 

The above results for // — > 777 can also be used to generate the helicity amplitudes 
for the bremsstrahlung processes 77 ff'j and e~7 £"77, if the crossing relations of 
Sect. ^]2| are applied. In particular, we get a set of helicity amplitudes for the so-called 
double Compton process 



e-(p,a) + 7(A;,A) e- {p' + j{k[, X[) + j{k'„ X',) 
by identifying / = e~, consistently replacing 



(5.17) 



P' 



//J- 



k 



3,A 



3,A, 
7 



3,A 
./ 

'i.A 



-k 



3,A' 



A 



-K 



a 



-a 



(5.1^ 



in all expressions, and adding the factor sgn((T'), where a' is the helicity of the outgoing 
e~ in ( |5.17|) . Finally, we mention that our helicity amplitudes have been numerically 
reproduced (up to phase factors) in the framework of the calculation presented in Ref. [|I^ 
by employing the spinor method of Ref. and by a third, completely different method. 
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Figure 4: Born diagrams for fi fi^ ^ f fli where graphs with the outgoing photon 
attached to the other charged fields are suppressed. 

5.3 The process — > f f") and related reactions 

As a final example, we consider the process 

f(q^r) + fiq',T')+^{k,\), (5.19) 

which will be of relevance at future muon colliders. The momenta p, p' and the corre- 
sponding Dirac spinors ^ for the incoming muons are defined in the same way 
as in the last section [see (|5.11|) and (|5.12 )]. The momentum matrices Q and Q' for the 



outgoing momenta q and q' are decomposed into the spinors pi^A and p[ ^, respectively, 

QAB = E P^,AP^,B, Q'ab = E P'iM,B- (5-20) 

1=1,2 1=1,2 

The respective Dirac spinors are generically given by 

^*' = (iMa), (5-21) 

with the actual insertions 




i.e. we have (^77) = {rj'C,') = ^f- The polarization spinors e* ^^l^) the outgoing photon 
are defined as in ( p.20|) . For f ^ p~ 16 diagrams contribute to the process in the Standard 
Model at tree level; they are schematically indicated in Fig. ^. The case f = n~ will not 
be considered here explicitly; the results for this particular channel follow from the ones 
for / 7^ /i~ by adding the amplitudes for the crossed process f p^ —>■ fp'^'j with negative 
sign. Since the interaction is mediated by the exchange of neutral particles only, the 
electromagnetic currents of the muon and the fermion / are conserved separately. This 
means that the gauge spinor for the outgoing photon can be chosen differently for initial 
and final-state radiation, which simplifies the calculation. In the following we generically 
denote the couplings of the neutral bosons to the fermions by 

-I- ^ -I- -^w ^ — ^ ^WJ 



g^f = -Qf, 9zf = -—Qf^ 9zf = -—Qf + ^ ' 

^± _ I riy_ , ^ ilwjriy 

fw' Mw 
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where Iwj = ±| is the weak isospin of the left-handed part of the fermion field /. In 
( ^.23| ) we follow the conventions of Ref. [jl2| for the Standard Model parameters and fields; 
in particular, x denotes the would-be Goldstone partner to the Z boson, and H is the 
physical Higgs boson. Since Standard Electroweak Theory does not conserve parity, the 
parity-induced relations between helicity amplitudes with opposite helicity configurations 
involve also an interchange of the chiral couplings: 



Mi-a, -a', -r, -r', -A) = sgn{aa'TT')Mia, a', r, r', A)^ 



(5.24) 



On the other hand, the process is CP-symmetric at the considered perturbative level, and 
the relations induced by CP symmetry read 



M(-a', -a, -r', -r, -A) = -sgnia a' tt')M (a, a\ r, r'. A)* 
The matrix elements 7V1 are decomposed according to the boson in the s-channel. 



A^((T,(T',r,r',A) = V2e^ 



E <Ha,a ,r,r 



(S), 



a, a' r, r'l 



(5.25) 



(5.26) 



LV=7,Z S=x,H 

leading to two generic functions A^^^ and A^^^ for the exchange of a vector field and a 
scalar field, respectively. Owing to (|5.24| ) or (|5.25|) , it is sufficient to give the results for 

A = +, 



4^)(a,a',r,r') 



(a,a',r,rO 



QpikP'Pk) 



QfikQ'Qk) 



2{p ■ k){p' ■ k)[{q + q'y- M^] 2{q ■ k){q' ■ k)[{p + p') 



vJ 



{q + q'Y - M2 



Qf 



^-^^^(9tmYm+9yf{kiY{H 
- [9ti,{kct>'Y{H + 9y,m*{i^'v)) 



QfikQ'Qk) 



Q^.{kP'Pk) 



4(g ■ k){q' ■k)[{p + p')^- Mi] A{p ■ k){p' ■ k)[{q + q'Y - M|] 

X [9UH')+9sM'^Y)(9tf{i'ri)+9sf{ivr 

1 



Qu9^,Ak^'Y{ki^Y / ^ , 



Qf9sf{kO*{kvr ,,,, 



+ 



'9t{w+9E,{<p'i^y 



2{p ■ k) 2{p' ■ k) 

1 1 

+ 



2(g ■ k) 2{q' ■ k) 



(5.27) 
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The quantities (kP'Pk) and (kQ'Qk) are defined in complete analogy to {kiP'Pki) in 
( CT ). The results for A^^ and Af^ are written down in the 't Hooft-Feynman gauge, 
i.e. we have = Mz- We mention, however, that we have reproduced the result for the 
complete matrix element A4 also in an arbitrary i?^-gauge for the photon and Z-boson 
fields, in which the individual contributions A^^'^ and A^^ look different. 

The above results have also been used to derive the helicity amplitudes for the process 



e-(p,a)+7(A:,A) e-(p', a') + e-(g, r) + e+(g', r') (5.28) 

in QED. The QED part of the amplitudes ( p. 261 ) is obtained by taking the contribution 
of V = 7 in (^.26|) only. Moreover, we identify = Qf = Qc = —1- The crossing 
proceeds analogously to the case described at the end of the previous section, i.e. one has 
to perform the replacements 

^ ^ kj^ ^ kj^^ k^ ^ k A ^ A, 

P'^^-P'^, ^'^A^-'^'^A^ o' ^ -o' (5.29) 



and to apply the sign factor sgn((T'), where a' is the helicity of the outgoing e~ in ( [5.28| ). 
Finally, we have to antisymmetrize all amplitudes with respect to the interchange of the 
two outgoing electrons. Moreover, the amplitudes have been numerically reproduced by 
an independent helicity method in the framework of the calculation discussed in Ref. . 



The polarized cross sections of e 7 — e e e+ contribute, in particular, to the left- 
right asymmetry Alr of polarized Compton scattering as background. The numerical 



agreement of the results presented in Ref. for this contribution to with the 



completely independent ones of Ref. represents an additional check of the calculation. 



6 Summary 

The Weyl-van-der-Waerden (WvdW) spinor technique for the calculation of helicity am- 
plitudes of massive and massless particles has been described in detail, providing all 
necessary ingredients for an implementation in computer algebra. This formalism leads 
to rather compact results for amplitudes, which can be directly used for numerical eval- 
uations, immediately after all spinor indices have been contracted to so-called spinor 
products. These contractions are technically similar to the usual ones for Lorentz indices. 
Since all mathematical objects, such as momenta, polarization vectors, and Dirac spinors, 
are expressed in terms of WvdW spinors, the spinor calculus often allows for further sim- 
plifications at the analytical level. Moreover, we have formulated how to exploit discrete 
symmetries for a reduction of the algebraic work or for providing additional checks, and 
how to avoid problems that are due to the appearance of unphysical poles in amplitudes. 

The use and the power of the described spinor technique have been demonstrated by 
the explicit calculation of the helicity amplitudes for the processes 77 //, // — > 777, 
\i~ [i^ — »• //7 with massive fermions, and of reactions that can be obtained from those by 
crossing symmetry. The results, which have been analytically simplified as far as possible, 
are very compact and well suited to numerical evaluations. 
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Appendix 



More details about discrete symmetries 



In this appendix we supplement the discussion of parity and CP symmetry in Sect. ^ 
by further details. 



(i) Parity: 

In order to explicitly derive the relations between helicity amplitudes that are con- 
nected by parity, we inspect the individual terms in an amplitude Ai after complex 
conjugation. Contractions between momentum matrices, which are just Minkowski 
inner products according to ( |2.13| ), are not changed at all, since the result is real. 
Contractions between a momentum k and a polarization vector e, or between two po- 
larization vectors, turn into the products for the corresponding polarization vectors 
with opposite hehcity, since (K^^)* = K^j^ = K^j^ and {s^^abT = ^\,ab = ^-\,ba 
(A = 0, ±1). Coupling factors iC for interactions between bosons simply turn into 
— iC*, and boson propagators receive a factor —1 owing to the explicit factor i in the 
numerator. The only non-trivial terms are those originating from Dirac matrices and 
fermionic wave functions. The complex conjugate amplitude remains unchanged if 
we multiply each fermionic vertex by a factor Fp = — 1 from the left and from the 
right, where the matrix Fp is defined as 











(A.l) 



Shifting then a factor Fp to the fermion propagators or fermionic wave functions 
that are attached on the right side of this vertex, and shifting a factor — Fp to the 
left side, we get the following replacements when going over from Ai to A4*: 



/2 



V/l/2 



i(|if + m/) 

;(±) 



(±) 
fc,2 



-Fp 
-Fp 
-Fp 

Fp m + rrif)] 



Fp = -ie^ (Cy/i/J 
^Fp = -ie^„, (C'vA/J" 
Fp = -ie (C^Ij^*uj^, 
Fp) = -i(^ + m^), 



Pp 
Tp 



k,l 
(±) 



k,2 



(±) 
fc,2 ' 

(±) 
fe,l ; 



(±) 



(±) 



-rp) 
-rp) 



(A.2) 



where the notation for the generic Feynman rules of Tables |l] and ^ is used. The 
relations between M.{ai] Xj)* and A^(— 0"^; — Aj) can be read off from ( [A.2| ) and the 
above considerations. The substitutions for the chiral couplings is explicitly given; in 
particular, there is no substitution in parity-conserving theories, in which Cyj^j^ = 

(pvfif'2} ^^'"^ There is a global change of sign due to vertices and propagators, since 
each of them introduces a factor —1, and since their total number is always odd in 
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tree-level amplitudes. Another factor —1 comes from our convention of extracting 
the factor i from the amplitude iAi. Finally, we encounter a factor sgn((jj) for each 
external fermion with helicity (Xj and a factor — sgn((Tj) for each external antifermion, 
leading us directly to relation (|4.3|) . 

CP symmetry: 

As in the case of parity, we derive the CP-induced relations between helicity ampli- 
tudes by considering the complex conjugate of a given amplitude A^. We assume 
an underlying model in which all couplings between bosons are CP-conserving, as it 
is the case in the Electroweak Standard Model. Then, for bosonic fields the situa- 
tion is simple, because the wave functions of particles and antiparticles are formally 
identical, and the bosonic couplings that are related by charge conjugation differ 
by at most the sign factors (see e.g. the Feynman rules in Ref. |]12|). Specifically, 
there is a factor (—1)" between each vertex factor iC and the expression iC for the 
vertex involving the corresponding charge-conjugated fields, where v is the number 
of vector bosons attached to the vertex. In order to get rid of the factor (—1)", we 
shift a factor —1 for each vector-boson leg of a vertex to the attached vector-boson 
propagators or external wave functions, resulting in an overall factor (—1)"^ for each 
diagram, where ny is the number of external vector bosons. After this procedure the 
bosonic propagators and vertices for the charge-conjugated fields consistently dif- 
fer by a factor —1 from the complex-conjugated counterparts involving the original 
fields. 

For fermionic fields we proceed in a way similar to the treatment of parity de- 
scribed above. In addition to taking the complex conjugate of a Dirac chain, we 
now also transpose each term and invert the order in the chain, since charge conjuga- 
tion reverses the directions of fermionic lines. Then we multiply fermionic vertices, 
propagators, and wave functions by the matrix 



CP 



1 
-1 



(A.3) 



with F^p 
ments 



-1 in the same way as done with Fp above. This leads to the replace- 



(±) 

k,2 



—Fcp 
—Fcp 
—Fcp 



-ieixCL 



Pcp = -ie|^ {Cyj^f,^* 
FcP = -ie^_A (C^/^j 



UO-r 



FcP = -ie {C-lJ 
i(^ + nif) Fcp W + ^f)]^ (-rcp) = + mj) 



t 



(±) 

fe,2 



-Fcp) 
-Fcp) 



(T) 
k,2 5 



k,2 



Fcp 
Fcp 



k,2 



■fl'g. (A.4) 



Note that we had to include a factor —1 for the vector-boson-fermion coupling 
in the first two lines, so as to be consistent with the above treatment of bosonic 
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couplings. Moreover, one should realize the correct change of sign in the momentum 
of the fermion propagator [see fourth line in ( [A.4|) ], corresponding to the inversion 
of the fermion line. In the case of CP-violation the necessary substitutions for the 
couplings can be read off from (|A.4|) ; if CP is conserved no substitution is necessary, 

because then CT.r r = (CTrr r ) and Clj r = (C^l r ) . The overall sign between 
Vht2 V V}ii2j •5/1/2 V •S/1/2/ ^ 

the amplitude M.{(Ji] and the one for the CP-related process M.{—ai] —Xj) is 
deduced as follows. At tree level, there is an overall factor —1 from the vertices 
and propagators and another —1 from the convention for Ai, as in the case of 
parity. From (|A.4|) one can see that incoming fermions and antifermions with helicity 
CTj introduce each a factor — sgn(cTj), while outgoing fermionic lines yield a factor 
sgn(crj). Taking into account the factor (— l)"^ derived above, we obtain relation 
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